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Abstract 

We prove a recent conjecture on the duaUty relation for correlation functions 
of the Potts model for boundary spins of a planar lattice. Specifically, we 
deduce the explicit expression for the duality of the n-site correlation func- 
tions, and establish sum rule identities in the form of the Mobius inversion 
of a partially ordered set. The strategy of the proof is by first formulating 
the problem for the more general chiral Potts model. The extension of our 
consideration to the many-component Potts models is also given. 
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I. INTRODUCTION 



The duality relation for the g-state Potts model is an identity |0 relating the par- 
tition functions of a Potts model on a planar lattice with that of its dual. Very recently, 
this duality consideration has been extended to the Potts correlation functions @J^. Specif- 
ically, it was established that certain duality relations exist for correlation functions of n 
Potts spins on the boundary of a planar lattice. Explicit expressions for the duality relation 
have been obtained for n = 2, 3 and n = 4 the expression for general n has also been 
conjectured 0. It has also been shown that there exist certain sum rule identities, but the 
explicit form of the identities was not given The purpose of this paper is the follow- 
ing: We first extend the consideration to the chiral Potts model and obtain its duality 
relation in a very general form. We then show that this formulation permits us to establish 
the conjecture of 0. In addition, it also leads to an expression of the sum rule identities in 
the form of the Mobius inversion of a partially ordered set. Furthermore, starting from a 
multi-component chiral Potts model ||^, we extend considerations to the multi-component 
Potts model, which includes the Ashkin- Teller model and its generalization This 
extends results reported elsewhere |Tl|] for the 2-component Potts model. 



II. THE CHIRAL POTTS MODEL 

Consider a g-state spin system on a two-dimensional lattice or, more generally, any 
planar graph C. Let Xi = 1, 2, ■ ■ ■ , g denote the spin state of the i-th site and —J{xi, Xj) the 
interaction energy between sites i and j, which can be edge- dependent. The interaction is 
chiral if 

J{x,x') ^ J{x\x). (1) 
The partition function of this spin system is 

E n Uix,,x,), (2) 

Xi = l <ij> 

where U{x,x') = exp[J{x,x')/kT] and the product is taken over all edges of C 

For our purposes we consider the chiral Potts model for which the Boltzmann factor 
is cyclic, namely, it satisfies 

U{x,x') = u{x — x'), (mod q). (3) 

We shall refer to the q x q matrix U as the interaction matrix. The chiral Potts model 
reduces to the standard Potts model upon taking J{x,x')/kT = K6{x,x'), where 6 is the 
Kronecker delta function which we define as 

6{x,x') = 1, X = x' (mod q) 

= 0, X x' (mod q). (4) 

This leads to 
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u{x) = l + {e^ -l)5{x,0). (5) 

Duality relation for the partition function: The partition function of the chiral Potts 
model possesses a duality relation 0. Let be the dual of C, and let Z^^\u*) be the 
partition function of the chiral Potts model on with the Boltzmann factoiQ 

u*{y) = ^Y^uj^'u{x), y = l,2,---,g (6) 

where uj = e^'^'^li. Then, the duality relation reads 

g-^/2Z[n(a:)]=g-^*/2z(^)K(y)], (7) 

where and A^* are, respectively, the numbers of sites of L and . They are related to 
the number of edges E by the Euler relation + A^* = £" + 2. This duality relation holds for 
arbitrary edge-dependent interactions. Note that y/qu*{y), y = 1, ■ ■ ■ , g, are the eigenvalues 
of the matrix U. In the case of the standard Potts model, (H) reduces to 



u*iy) = ^ 



e^-l + g<5(y,0) 



(8) 



III. DUALITY RELATION FOR THE CORRELATION FUNCTION 

Number n sites on the boundary of C by integers 1, 2, ■ ■ - n clockwise as shown in Fig. 1. 
The probability that the n sites are in respective spin states xi,X2, ■ ■ ■ , x„ is 

Pn{Xi, X2, ■ ■ ■ , Xn) = Z^^^^...^^/ Z (u) , (9) 

where Zx-^^x2---xn is the parim/ partition function of the chiral Potts model with the n spins in 
fixed given states. Construct an auxiliary lattice £aux from C by connecting all n boundary 
sites to an extra site e in the infinite face of C as shown. These n connecting lines then 
divide the infinite face of C into n parts, which we also number by 1,2, - ■ ■ ,n as shown in 
Fig. 1. Let Uje, j = 1,2, ■■ - n, be the Boltzmann factor of the line connecting sites j to 
e. Its dual edge then connects sites (residing in faces of C) j and j + 1 on with the 
Boltzmann factor 

- Vj) = ^ E M,e(x), J = l,---,n (10) 

V ^ x=l 

where yj denotes the spin state of the dual site j. Here we have adopted the convention of 
in orienting the edges. It is also convenient to consider the n-spin correlation 0| 



^The definition of the dual Boltzmann factor adopted here differs from that used in by a factor 
y^. The present definition is consistent to the notion of (u*)* = u. 
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(11) 



a quantity which vanishes identically if the n spins are completely uncorrelated. 

We apply the duality relation to >Caux- Let Z*^y^,,,y^ be the partial partition function 
of the Potts model on with sites 1, 2, ■ ■ ■ , n in respective states yi, y2,---, Un- Then, 
explicitly, (0) reads 



^ {N+l)/2 ^Uie{Xie)u2e{x2e)---Une{Xne)ZxiX2-x,, 

XlX2---Xn Xe 

= q-iN*+n-l)/2 ^ ^,(y2l)^32(y32)---<(2/ln)^;,,...,„, (12) 



yiy2---yn 



where Xie = Xi — Xe, 1/21 = 1/2 — Vi, etc., and we have used the fact that £aux has + 1 sites 
and has N* + n - 1 sites.0 

Due to the cyclic nature of the matrix elements, we can replace the arguments Xie, X2e, ■ ■ ■ 
on the LHS of (0) by Xi,X2, ■ ■ ■■ The summation over x^ can then be carried out leading 
to an overall factor q. Next we eliminate u* by using (0). Since all u factors in (|12]) are 
independent, as a consequence we can equate the coefficients of Uie{xi)u2e{x2) ■ ■ -Uneixn) in 
(|1^). This leads to the duality identity 

^xix2-x„ = C'„(g) M{xi,X2,---,Xn\yi,y2r--,yn)Z*^y^...y^, (13) 

yiy2---yn 

where 

Cniq) = (14) 

and 

M{xi,X2, ■■■,xn\yi,y2,---,yn)= (15) 



can be interpreted as the elements of a x g"' matrix M. Eq. (13) is the duality relation 
for the chiral Potts model, and is a very general expression. The relation assumes a reduced 
form for the standard Potts model after taking into account the degeneracy of states, and 
the crux of matter is, of course, to deduce the reduced relation. This is the subject matter 
of the next section. 

Introducing the notation 

Py.y.-y. = ^^..-.J^n...! = ^Z;^y.-yJZ^''^ (16) 
and combining with (|^) and (P), we can rewrite (p!3|) as 

P„(a;i,a;2,- ■■,x„) = J2 (xi, X2, ■ ■ ■ , x^-i, | ?/2, ■ ■ ■ , l/n)Pyi2/2-2/n- (17) 

3/13/2 ■■■yn 

An immediate consequence of ([T7|) is 



^Note that the partial partition function Z*_^y^,„y^^ in (12) differs from that in yj5[ by a factor. 
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Tn ^ Pyiy2---yn 1) (18) 

" yiy2--yn 

a result which has been reported previously |§|,|12| . 

To obtain the inverse of (|T3D , we note that, by deleting the {1,2}, {2,3}, ■ ■ ■ edges in 
^aux5 arrive at a lattice C* in which 1, 2, ■ ■ ■ of are its boundary sites. An example of 
£ and the resulting C* is shown in Fig. 2. It is clear that the roles of C and C* are reciprocal, 
so we can apply (|T3|) to C* to obtain an inverse transform. Now C* has N = N* + n — 1 
sites and its dual has N* = N — n + 1 sites. As a result, we obtain 

^yiy2-y« = ^niq) J2 ?/2, " " " , 2/n | Xi, Xa, ■ ■ ■ , a;„_i)Z^^^2...^„, (19) 

xiX2---x„ 

where 

C„(g) = ^-n+(7V-iV*)/2 ^ ^-l+(^*-7V)/2_ (20) 

We have c„(g)c„(g) = and note that the sequence {xn,xi,X2, ■ ■ ■ ,Xn~i} of the in- 

dices on the RHS is related to {xi, X2, ■ ■ ■ , Xn-i, x„} by a simple cychc permutation. Upon 
combining (|19D with (|13D we obtain the identity 

M^(Xi, X2,---,Xn\ X[, X2, ■ ■ ■ , x'J = q^'^^6{xi, X2)S{X2, Xg) ■ ■ ■ x[) , (21) 

which we refer to as a reciprocal inversion relation. 



IV. THE STANDARD POTTS MODEL 

For the standard Potts model the duality relation (0) assumes a reduced form which 
has been given previously for n = 2, 3 0] and for n = 4 0. The form for general n has also 
been conjectured 0. It has also been shown that the correlation functions satisfy certain 
sum-rule identities 0. Here, we reformulate the standard Potts model as an instance of the 
chiral Potts model and, as we shall see, this leads us to establish the conjecture as well as 
deduce the general expression for the sum rule identities. For this reason, it is instructive 
to first demonstrate how ([T3| ) reduces to the known results for n = 2, 3. 

The standard Potts model is characterized by the fact that one needs only to keep track 
of spin states that are the same. That is to say, one needs to keep track of the partition 
of the n sites into blocks, such that sites in one block are in the same state. This leads to 
writing the partial partition functions in an expansion of a partially ordered set. Along this 
line, we writes for n = 2, 3, 

Du + Du6{xi,X2) 
Dl, + Dl,6{yi,y2) 

Dl23 + Du35{Xi,X2) + Di2lS{Xi,Xs) + Di22S{x2,X3) + DiuS{Xi, X2, X3) 

D*i23 + Dl,s5{yi, 2/2) + Dl2^5{yi, y^) + 1^^22^(2/2, ys) + ^ni5(2/i, 1/2, Vs), (22) 

where 6{xi,X2,X2,) = 6{xi,X2)6{x2,X3). Here, for n = 3, for example, the partition of the 
3-set of integers {1,2,3} includes the 5 elements {{1}, {2}, {3}}, {{12}, {3}}, {{1},{23}}, 
{{13}, {2}}, {{123}} denoted by subscripts {123}, {113}, {122}, {121}, {111}, respectively. 



^XlX2 

yiyz 



^XlX2X3 



yiy2y3 
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For n = 2, we substitute the first two lines of ( ]22|) into (|T3]). Making use of the identity 



y=l 



(23) 



one arrives at 



Du + DnS{x,,X2) = C2{q) E oo^''^'^^''^' 



C2{q) 



q'Dl^6{xi,X2) + qDl^ 



(24) 



Since this equation holds for arbitrary xi and X2, the coefficients of corresponding delta 
functions must be equal, and this leads to the identities 



Di2 = C2{q) qDl 



ll5 



Dn = C2iq) q'Dl,, 



(25) 



which can further be converted to relate Zx^x2 Z*_^ as follows. First, the second line of 



gives 



^12 ~ -'^125 



Zl, = Dl, + Dl,, 



from which one obtains the inverse 



^12 — ^125 



D 



11 



^11 



'12- 



(26) 



(27) 



Substituting (^) into (^51) and combining (H), (|^) and the definition (|T6D, one is led to the 
duality relation for the 2-point correlation function 



P2{Xl,X2) = [Di2 + Du6{Xi,X2)]/Z 

= — 1 + (^q5{xi, X2) - 1 )pi2 



This is the duality relation ([T7|) for n = 2. Finally, using (p2D, (p5D, and 
the row vector Z2 = (^11, Z12), = {Z*^, ZI2), we find 



Z2 = gC2(g)T2(g) ■ z^. 



(28) 

), and introducing 
(29) 



where 



T2(g) 



1 q-l 
1 -1 



(30) 



satisfying [T2(q')]^ = qi-21 with I2 being the 2x2 identity matrix. 
Similarly for n = ?> one obtains in place of (0) 

^123 + Dii^5{Xi,X2) + L'i2i(5(Xi,X3) + £'l225 (a;2, ^^s) + £'lll5(Xi,X2,X3) 



CM E 



a:31?/l+3;i2J/2+2;23?/3 



^123 + ^113^(?/1' ?/2) + ^t21^(?/l> ?/3) 



?/l?/2?/3 

+^i22'^(l/2, 2/3) + ^ni'^(2/i> 2/2, 2/3) 



(31) 
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leading to the identities 



^123 = CM qDul, Dm = OM q'Dl,, 

Du3 = Csq'Dl,„ = C3(g) g'/^t22> Du2 = Cs{q)q'D 



2 n* 
113- 



(32) 



Again, we can rewrite (|32D in a form relating Zxix2X3 to Z*^y^y^. Now the inverse of the last 
line of (521) is 



-'-^123 


— ^123 






-'-^113 


— ^113 ~ 


ry-^ 7~j* fy-if. ry-if- 7^* 

^123' -^121 ~ ^121 ~ ^123' -^^122 ~ ^122 


~ ^123' 


^111 


— -^111 ~ 


(^113 + ^121 + ^122) + 2^r23- 




Substituting (^) 


into (H) 


and using the third line of (p2D for Zx^^^xz-, 


we obtain 



(33) 



1 - (Pl22 + P12I + P113) + 2pi23 + g(Pl21 - Pl23)S{Xi, X2) 



+g(Pl22 - ^123)^(3^1, 3^3) + q{Pll3 - P123)S{X2, X3) + q'^Pl23S{Xi, X2, X3) 



(34) 



This is the duality relation (|T^ for n = 3 reported in 

Explicitly, introducing the row vector Z3 = {Zm, Zu^, Z121, Z122, Z123) and similarly 
defined Zg, we find 



Z3 = qC3{q)T3{q) ■ z 



3; 



(35) 



where 



T3(g) 



satisfying 



/I q-1 


1 


- 1 


q 


- 1 


{q 




i)(?-2: 


\ 


1 -1 


q 


- 1 




-1 




-(g-2) 




1 -1 




-1 


q 


- 1 




-(g-2) 




1 q-1 




-1 




-1 




-(g-2) 




Vl -1 




-1 




-1 






2 


/ 






/I 






























1 









[T3{q)r = 







1 

























1 
















Vo 











l) 







(36) 



(37) 



^Note, however, the numbering of {2/1,2/2)2/3} as {s',s",,s} in resulting in a more symmetric 
appearance of (E 
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V. THE GENERAL ANALYSIS 



For general n, we write Z^^-.-xn as a sum in the form of (0), namely, 

Zx,...x„ = Y.Dx5{X), (38) 



X 



where the summation is over all partitions X of n integers {1,2, ■ ■ ■ ,n}, and 6{X) is a 
product of delta functions associated with blocks in X. For example, we have S{X) = 
6{xi,X2)S{x3,X4) for X = {1133}, and S{X) = S{xi,X2,X3) = 6{xi,X2)S{x2,X3) for X = 
{1114}. If a block of X contains m > 3 x's as in S{xi,X2,X3) above, we write S{X) as a 
product of m — 1 2-point delta functions 6{xj,Xk) with the x's arranged in, say, a clockwise 
sequence, so the writing of S{X) as a product of 2-point delta functions is unique. The 
number of terms in (^) is 6„, and it has been shown in that bn is generated by 

oo 

^M7^! = exp(e*-1). (39) 

n=0 

Without the loss of generality we shall assume q > n so that all Zx are physically realized.^ 
The partition sums (E2|) and (pH]) can be regarded as a transformation between the Z's 



and D's. The transformation is that of a partially ordered set and it is known |]13l that 
the inverse is given by the Mobius inversion. Let X, X', Y, etc. denote specific partitions 
of an n-set. We write X' ^ X if the block structure of X' is contained in X, namely, X' 



is a refinement of X For example, {1134} and {1234} are contained in {1133}, while 



{1224} is not. Then, an immediate consequence of (^) is ||13| 



?x = E Dx'. (40) 

X'^X 



In addition, the Mobius inversion of (^) is 

Dx= Yl KX',X)Zx' (41) 



X 

with 



MX',X) = (-1)I^'H^I n K-l)!> (42) 



blocks eX 



where |X| and \X'\ are, respectively, the number of blocks in X and X', and rib the number 
of blocks of X' that are contained in a block of X. In the examples above, for instance, 
X = {1133} has two blocks, and we have = {1,2} for X' = {1134}, and rih = {2,2} for 
X' = {1234}, etc. This leads to 



^Otherwise the situation q < \X\ < n can occur, where \X\ is the number of blocks of the 
partition X, and Zx has no physical meaning. But Zx, which can be computed for q > \X\, is a 
polynomial in q. We shall take Zx to be the same polynomial for all q including q < \X\. With 
this understanding in mind all our results, including the identities derived below, hold for all q, n. 
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1^1133 = (-l)'-'(0!)^Zn33 + (-l)='-'(0!)(l!)Zn34 + (-l)='-^(l!)(0!)Zi233 + {-lf-\VfZ,2M 

= ^1133 ~ ^1134 — ^1233 + ^1234- 

(43) 

Other examples are 

DlllA = Ziiu — ^1134 — Zi2lA — ^1224 + 2^1234 

-D1212 = ^1212 — ^1214 — ^1232 + ^1234- (44) 

Substitute the partition sum ( ^Sf ) and a similar expression for Z* .. into (|1^), we obtain 



J2Dx6{X) = C^{q)J2 E M{x,,---x^\y,,---yn)D*y6{Y). (45) 

X Y yi- -yn 

It turns out that the analysis is best done graphically. Represent the matrix element 
M{xi, ■ ■ ■ ,Xn\yi, ■ ■ ■ ,yn) in (0) by a graph G„ shown in Fig. 3, where each node (open 
circle) denotes a ?/-summation, and each arc is assigned a label x such that an outgoing 
arrow (from a node labelled y) carries a factor uj~^y and an incoming arrow a factor u^^. To 
each X we construct a "connectivity" Tx according to the following prescription: Connect 
the mid-points of the arcs belonging to each block of X to a common point exterior to Gn- 
A connectivity is planar if the connecting lines do not intersect; otherwise the connectiv- 
ity is non-planar. Examples of connectivities for G4 are shown in Fig. 4. We note that 
the construction of the connectivities is the same as that introduced in [^]. We have two 
possibilities. 

i) Tx is planar. In this case Tx divides the region exterior to Gn into faces. Regarding 
all nodes in one face as belonging to a block of another partition Y, and denote this mapping 
hj X Y. It is clear that Fy is also planar and that Y ^ X. Therefore the mapping 
X y is one-one. One has alsoQ 

|X| + |F|=n + l. (46) 



It has been shown that the number c„ of planar Fx, or simply planar X, is generated 

by 

00 

E c„r = (1 - VT^)/2t. (47) 

n=0 



With these understandings, we now carry out the y-summations on the RHS of (^51). 
From (^) each y-summation yields a factor q 6{x, x'), where x and x' are the labellings of 
the incoming and outgoing arrows. One can then regard the labellings x as a flow which is 
conserved at each node. Further, the effect of the delta function 6{Y) in (^) is to collapse the 
y-summations of each block of Y into a single summation, resulting in a factor q 6{xia, Xout) 
where is the sum of the incoming flows and Xout the sum of outgoing flows. This effect 



^For X = {{a;i},{x2},---,{x„}}, |X| = n, we have Y = {{yi, ^2, • • • , 2/n}}, \Y\ = 1. It is also 
clear that |y| increases by 1 when |X| decreases by 1. This establishes 
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can be conveniently visualized by contracting the nodes in question into a single one, and 
requiring that the flow is conserved at the contracted node. For example, as shown in Fig. 
5, the partition Y = {123216} of Gq is contracted into a diagram consisting of a sequence of 
3 G2 "cactus leaves" attached together at two common nodes. Generally, the delta function 
S{Y) contracts On into a diagram consisting of cactus leaves Gm, m < n, attached together 
at common nodes. The planar nature of Y then ensures that the resulting diagram is a 
cactus tree. One next carries out the y-summations one-by-one starting from the nodes of 
the outermost leaves of the tree. The tree structure now ensures that this always leads to a 
product of 2-point delta functions. In the example shown in Fig. 5, for instance, this leads 
to 

5(X5, XG)S{Xi + X6, Xi + X5)5(Xi + X3, X2 + X4)5{x2, X3) = 5{x5, XG)5{Xi, Xi)5{x2, X3) . (48) 

But this is precisely the factor 6{X) appearing in DxS{X) on the LHS of (^), where 
X = {122155} Y. It is readily verified that, generally for each planar X ^ Y, to the 
term DxS{X) on the LHS of (^), the ?/-summations of the Y term on the RHS yields a 
factor Dygl^l5(X). Now, the indices {xi} in a given S(X) are arbitrary. It follows that the 
coefficients of the two terms must be equal, and this leads to the identity 

Dx = Cn{q) q^^^D*y, planar X ^ F, (49) 



an equation which generalizes (^) and 

ii) Tx is non-planar. We now consider the remaining 6„ — c„ terms in (^5]). Now, each 
remaining terms on the LHS is of the form Dx8{X), where 6{X) is a product of 2-point 
delta functions ^(xj, Xfc). For the remaining terms on the RHS we again carry out the y- 
summations for each Y . However, due to the fact that Y is now non-planar, the process of 
contraction as described in the above leads to diagrams in which cactus leaves form rings, or 
circuits. It follows that the y-summations generates a product involving some 4-point delta 
functions of the type 5{xi + Xj, Xk + Xi), which cannot be reduced into products of 2-point 
delta functions as in (^8]) . In other words, the delta functions 5{X) occurring on the LHS 
do not appear on the RHS, and vice versa. Now, the the indices {xj} in a given 6{X) are 



arbitrary. It follows that, for the equality (45) to hold for all {xj}, each of these remaining 



coefficients must vanish individually. Namely, we must have 

Dx = D*Y = Q, non - planar X, Y. (50) 
Finally introduce the partially ordered partition sums 

P„(xi,X2, ■ ■ ■ ,x„) = ^Ax6{X) 

X 

Py.y.-yn=T.BY6{Y), (51) 



Y 



where Ax = Dx/Z, Bx = Zp/Z 11...1- Combining (|^), (|l^) and using (|46|), the identities 

(^) and ( 1^01) now establish the conjecture of [Q, which we now state as a 

Theorem: 

Ax = q'^^^By, for planar X ^ Y, 

= 0, otherwise. (52) 
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Explicitly, the coefficients By are given in terms of Pyiy2---yn by the Mobius inversion (^TJ). 
Thus, the second line in ( ^2|) leads to identities relating partial partition functions Zx- For 
example, from the second line of (0) and A1212 = -D1212 = 0, we obtain 

Z1212 — 2'i214 — ^1232 + ^1234 = 0. (53) 

This is the n = 4 identity reported in |Q. Generally, one obtains an identity Ax = for 
each non-planar X. Using these equations one can express all non-planar Zx in terms of 
planar ones, a fact reached in through the use of high-temperature expansions. 

One can also write down the transformation relating Zx and Zy- For general n and 
planar X one finds by combining (|40|), (^) and (^T]), 



Zx = Cn{q) E E f^{Y",Y')z;.., x'^r. (54) 



\y' 

X'<X Y"-<Y' 



The transformation relating Pn{xi, ■ ■ ■ ,Xn) to Py^-y„ can be written down similarly with 
the factor C„(g) replaced by More explicitly, define a row vectors z„ and z* whose 

elements are the c„ partial partition functions corresponding to planar X, we find 

z„ = gC„(g)T„(g) ■ z^, (55) 

where T„(g) is a c„ x c„ matrix satisfying the identity 

[T„(g)]2(X, X') = g"-i5(xi, 4)5(X2, 4) " " " ^{^n, x[). (56) 

Note that the coefficient on the RHS differs from that in (0). Expressions of T„(g) for 



n = 2, 3 have been given in ( |5(]| ) and (^). The expression for n = 4 can be deduced from 
results reported in 0. For n = 5 there are 42 planar X and we define the row vector 



Z5 — ( -Z^lllll; (-2^11115; -2^11141; -2^11311) ^12111! ^12222); (^11144) ^11331) ^12211; ^11333 


(^11145) 


^11341; 


^12311) ^12333) ^12225) 


) (^11335; ^12241; ^11344? ^12331) 


^12244); 


(■Z^11315j 


2'l2141, 


2^12311; 2'i2115, Z122A2) 


J (-2^11343; 2^12321, -2^12144, -2^12332, 


■2^12215); 


(^12125) 


^12142) 


^12312) ^12313) ^12323) 


) (^113457 ^12341; -2^12344! ^12335) 


^12245)5 


(^12145) 


^12342; 


^12315? ^12343) ^12325) 


5 ^12345 ) • 





(57) 



Then after some algebras we deduce from (Q) the explicit transformation 

^11111 = {1 + gi(i, 1, 1, 1, 1) + gi(i, 1, 1, 1, 1) + qMi, 1, 1, 1, 1) + qMi, i, i, i, i) 

+g?(l, 1, 1, 1, 1) + ql{l, 1, 1, 1, 1) + gis(l, 1, 1, 1, 1) + qlq2{l, 1, 1, 1, 1) + gigaO, 

^11115 = {1 + (-1, qi,quqi, -i) + (-1, qi, qi, -1, qi) - ^2(1, -gi, -gi, 1, 1) 

-^2(1, -qi, 1, -qi, 1) - gi(i, -qi, 1, 1, 1) - gi(i, -gi, 1, 1, 1) + s(i, -gi, 1, 1, 1) 

-gig2(l, 1,1,1,1) -g2t}, 
^11144 = {1 + (gi, -1, gi, gi, -1) - (1, 1, -gi, 1,1)- g2(i, 1, -gi, 1, 1) - g2(i, 1, 1, 1, 1) 

+(-gi, -gi, -gi,gi, 1) + (i,-gi,-gi,i,g?) - ip,s,s,p,p) 

+g2(-gi, 1, gig2, 1, -gi) - g2M}, 
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2'iii45 = {1 - (1, 1, -qi, -qi, 1) -(1,1, -qi, i, i) + (2, -g2, gig2, -g2, 2) - (-2, gs, g2, g2, 92) 

-(gi, gi, gi, gi, -1) - (-1, gi, gi, -1, gi) + (/i, -s, -s, h, h) 

+(2gi, g2, -gig2, ^2, 2gi) + 2(t + g)}, 
^11335 = {1 - (1, -gi, 1, -qi, 1) - (1, 1, 1, 1, 1) - g2(l, 1, 1, 1, 1) + (2, -g2, 2, -gs, 2) 

+ (1, 1, -qi, -qi) + {-qi, l, l, l) + (^, -P, h, h, h) - (gigs, g2, g2, gig2, 2) + 2m}, 
^11315 = {1 - (1, 1, 1, -gi, 1) - (1, -gi, 1, 1, -gi) - (-2, g2, g2, g2, 92) + (-92, -g2, 2, gig2, 2) 

+ (1, 1, 1) + (1, -gi, 1, -91, 1) + (e, -s, /i, -s, e) + (2gi, ga, 2gi, 2gi, gs) + v}, 
^11343 = {1 - (1, 1, 1, -gi, 1) - (-gi, 1, 1, 1, 1) - (<?2, -2, g2, g2, -2) - (-2, -2, g2, g2, g2) 

+ (1, -gi, 1, -gi, 1) + (1, 1, 1, 1, 1) + (e, -p, h, e) - (gig2, 2, -ga, -g2, 2) + i; + 3g}, 
^11345 = {1 - (1, 1, 1, -gi, 1) - (1, 1, 1, 1, 1) + (2, 2, -ga, -g2, 2) + (2, 2, 2, -ga, 2) 

+ (1, -gi, 1, -gi, 1) + (1, 1, -gi, 1, 1) + (-5, h, h, h, -5) + (2gi, -2, ga, g2, -2) + 14 - 5g}, 
^12145 = {1 - (1, 1, 1, 1, 1) - (1, 1, -qi, 1, 1) + (2, 2, -ga, 2, 2) + (2, -q^, 2, 2, -ga) 

+ (1, 1, 1, 1, 1) + (1, -gi, 1, 1, -gi) - (g + 5, 5, 5, 5, ~h) - (2, 2, -ga, 2, 2) + 14 - 2g}, 
^12345 = {1 - (1, 1, 1, 1, 1) - (1, 1, 1, 1, 1) + (2, 2, 2, 2, 2) + (2, 2, 2, 2, 2) 

+ (1, 1, 1, 1, 1) + (1, 1, 1, 1, 1) - (5, 5, 5, 5, 5) - (2, 2, 2, 2, 2) + 14}, (58) 

where g^ = q—m, e = g — 5, h = 2g— 5, s = g^ — 5g+5, p = g^— 4g+5, t = g^ — 7g+7, u = q^ — 
5g+7, v = g^ — lOg+14. Here, one reads off elements of T5(g) from (^) directly. For example, 
the second row of T5(g) has elements {1, (—1, gi, gi, gi, —1), (—1, ■■■),■■■, —q2t}- Elements 
not shown are given by cyclic permutations of the partition indices within each set of paren- 
theses. For example, the third row of T5(g) is {1, (—1, —1, gi, gi, gi), (gi, ■■■),■■■, —g2t}, 
where elements within each set of parentheses are obtained by cyclically permuting those in 
the second row. 

VI. THE MULTI-COMPONENT POTTS AND CHIRAL POTTS MODELS 

The chiral Potts model can be generalized to more-than-one components. An m- 
component chiral Potts model is a 

Q = N^N2---N^ (59) 

state spin model, where N^^C. = 1,2, ■■■,m, are positive integers. The interaction matrix 
\J is a Q X Q matrix in the form of a direct product of m cyclic matrices, with elements 
indexed hj U{xi, ■ ■ ■ ,Xm\x[, ■ ■ ■ , x'^), xi, x'^ = 1,2, ■ ■ ■ , Ni, i = 1, ■ ■ ■ , m, satisfying the cyclic 
property 

, Xm I ■ ■ ■ , X^^ u(^Xi Xi, ■ ■ ■ , Xjn -^m)' 

(60) 

The m = Ni = N2 = 2 version is known as the Ashkin- Teller model 0. 

It is convenient to introduce vectors x = (xi, ■ ■ ■ ,Xm) and y = {yi, - ■ ■ ,ym) and their 
scalar product 

X- y = xi?/i/A^i H \-Xmym/Nm, X i, y^ = I , ■ ■ ■ , N i. (61) 
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Then the dual model has the Boltzmann factor 



Ni N„ 



^ ■ ■ ■ ^ e"""" M(x) 



VQ xi=i 

and the duality relation (0) now reads 

Q-^/2z[«(x)] = g-^*/2z(^)K(y)]. 



(62) 



(63) 



Here, again, y/Qu*{'y) are the eigenvalues of the interaction matrix U. 

Proceeding in a similar fashion as before, we arrive at the following duality relation in 
place of (p!3D, 



CniQ) J2 ^^(xi'X2, ■ ■ ■ ,x„ I yi,y2, ■ ■ ■,y„)Z*^ 



y2-yr7 



(64) 



yiy2---yn 

Here, M is a x matrix with elements 



M(xi, X2, ■ ■ ■ , x„ I yi, y2, ■ ■ ■ , yn) = exp 



27ri(x„i ■ yi + xi2 ■ y2 H h X(„,_i)„ • y„) 



satisfying the reciprocal inversion relation 

M^(xi,X2, ■ ■ - .Xn I xi,X2, ■ ■ ■ ,X^) = Q"+^(5(xi , X2)(5(x2 , X3) ■ ■ ■5(x„,Xi). 

Defining the n-point correlation functions 

P„(xi, X2, • ■ ■ , X„) = ZxiX2..-x„/^(m), 



(65) 



(66) 



(67) 



then ( p^ ) becomes 

P„(Xi,X2,---,X„) =Q-("+l) ^ M(Xi,X2,---,X„ |yi,y2,---,yn)Pyiy2-yn- 

yiy2---yn 

This is the most general correlation duality relation for the m-component chiral Potts model. 



(68) 



An immediate consequence of ( |68D is, in analogous to (11) and (0), the following expression 
for the n-point correlation 



r„ = Q'^P„(x,x,---,x)-l = g-i Pyiy2-y„ - 1- 



(69) 



yiy2---yn 



The m-component Potts model: The m-component chiral Potts model reduces to an m- 
component (standard) Potts model when the interaction ( |60|) assumes the form (^) for each 
component, namely, 



UyXi , ■ ■ ■ , Xm J 



n 



l + (e^^-l)5;v,(a;£,0) 



l,2,---,iV,, 



(70) 



where 
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SNi{x,x') = 1, X = x' (mod Ni), 

= 0, xj^x' (mod Ne). (71) 



The m = 2 version is known as the {Ni, N2) model . 



Now each Xj in (|68D is an m-component vector. Let its components be Xj£ = 1, ■ ■ ■ , A'^, i = 
1, ■ ■ ■ ,n; i = 1, ■ ■ ■ , m, and let Xi denote the partition of the n integers 1, 2, ■ ■ ■ , n dictated 
by {xii, i = 1, - ■ ■ , n}. Then, in analogous to (|5T|), we write the partially ordered mappings 



P„(xi,X2,- ■■,X„) =J2^X^X2-Xm^NA^l)SN2{X2) ■■■SNmi^m), 

Py.y.-y. = Bny2- YjNAYi)SN2iY2) ■ --SM^iY^). (72) 



With these notations, we have the following corollary of our Theorem: 
Corollary: 

For the m-component Potts model whose correlation functions are (|72D, we have 



Ax,X2-Xm = I[Ne By,y2-y^, for planar ^ Ye, 

£=1 

= 0, otherwise. (73) 



Here, explicitly, by using 

ZxiX2-Xm/Z = X! ^X[X'2-X^ 



X'g-<Xi 



BnY2-Y^ = E Pnn -y^ E -"(^/^ ^e), (74) 



in analogous to (0) one can reduce the correlation duality relation (0) to 

ZxiX2--X,^ = Cn{Q) E E ^Yl'Y^'-Y^ 
X[-<XiY'^<Y'^ 
m 

xn^VrWA^^/), X'.^Yl. (75) 
e=i 

This is the desired duality relation for the m-component Potts model. Now there are (c^)™' 
partial partition functions Z and Z*. Considering Z and Z* as tensors and of rank 
m, one can rewrite (^) more compactly as 



Zm — QCn{Q) 



T„(iVi) ® ■ ■ ■ ® T„(iV„ 



■ Z:;, (76) 



where the tensor products are over m Cn ^ Cn matrices T„(g) defined in (0). The last 
expression generalizes the m = 2 results for the (A^i, N2) model reported in |[Tl[] . 
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VII. SUMMARY AND ACKNOWLEDGEMENT 



We have considered the chiral Potts model and obtained the duahty relation for its 
correlation functions of n sites on the boundary of a planar lattice. The result is given by 
(p!3|) and ([T7|) . By specializing this result to the standard Potts model, we establish a recent 
conjecture of on the correlation duality which we now state in (|52D as a Theorem. The 
explicit duality relation relating the partial partition functions is given in (|5^) , with the 
n = 2, 3, 5 expressions explicitly worked out. The formulation is next extended to the multi- 
component chiral Potts model, leading to the correlation duality relations (|6^) and (|68|). 
Again, specializing the results to the m-component Potts model, we obtain the Corollary 
([73|) and the correlation duality relations (|75D and ([76|). 

Work has been supported in part by NSF Grant DMR-96 14170. 
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Figure captions 

Fig. 1. A planar graph L and n sites i, j, ... m, £ on tlie boundary. 
Fig. 2. Reciprocal graplis L (solid lines) and C* (broken lines). 
Fig. 3. Graphical representation G„ of the matrix element (|65|) . 

Fig. 4. Examples of connectivities on G4. (a) X = {1133}. (b) X = {1114}. (c) X = 
{1212}. 

Fig. 5. The contraction of the the Y = {123216} graph of Gq into a cactus tree. 



17 



Fig. 



I 1 



\ 



— X — 



— X 



X — 



- X 



\ 
X 



y4 



Fig. 2 




Fig. 3 



(a) 



(b) 



(c) 



Fig. 4 




(a) (b) 



Fig. 5 



